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ABSTRACT 


Context. Superbubbles and supershells are the channels for transferring mass and energy from the Galactic disk to the halo. Magnetic 
fields are believed to play a vital role in their evolution. 

Aims. We study the radio continuum and polarized emission properties of the W4 superbubble to determine its magnetic field strength. 
Methods. New sensitive radio continuum observations were made at 46 cm, 411 cm, and 421 cm. The total intensity measurements 
were used to derive the radio spectrum of the W4 superbubble. The linear polarization data were analysed to determine the magnetic 
field properties within the bubble shells. 

Results. The observations show a multi-shell structure of the W4 superbubble. A flat radio continuum spectrum that stems from 
optically thin thermal emission is derived from 1.4 GHz to 4.8 GHz. By fitting a passive Faraday screen model and considering the 
filling factor f Ue , we obtain the thermal electron density n e = 1.0/ y[f^ e (±5%) cm -3 and the strength of the line-of-sight component of 
the magnetic field 5// = -5.0/ yjf^ e (+10%) pG (i.e. pointing away from us) within the western shell of the W4 superbubble. When 
the known tilted geometry of the W4 superbubble is considered, the total magnetic field B tot in its western shell is greater than 12 pG. 
The electron density and the magnetic field are lower and weaker in the high-latitude parts of the superbubble. The rotation measure 
is found to be positive in the eastern shell but negative in the western shell of the W4 superbubble, which is consistent with the case 
that the magnetic field in the Perseus arm is lifted up from the plane towards high latitudes. 

Conclusions. The magnetic field strength and the electron density we derived for the W4 superbubble are important parameters for 
evolution models of superbubbles breaking out of the Galactic plane. 
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1. Introduction 


A large number of superbubbles and supershells have been 
found and discussed with i n our own Galax y (e.g. HeBei 1979; 
Hul Il98l1: iHeilesI 1 1984 iKoo et al.l Il992l: [ Macieie wski et al.l 


lUvamker et al 


20011: iMcClure-Griffiths et al.1 120001 

120071) and in external galaxies (e.g. 


Pidopryhora et al 


Meaburn ll980l:lGraham & Lawrielll982l:lBrinks & Bajajal fl 986; 


Deul & den Hartogll 19901: iKim et alJl 1998b . They are created ei 

ther by energetic stellar winds from OB star associations, by 
multiple supernova explosions, or by a combination. Besides 
their identification by shell-like HI emission, superbubblesand 
supershells were also t raced by soft X-ray (e.g.lCash et al.lll9~80h 
and H a emission (e.g. [Reynolds & Ogden![19791) . Superbubbles 
power ed by sufficient energy wi ll break out of the Galactic 
plane dTomis aka & Ik euchilTl986l: iMac Low & McCravl Il988l: 
Mac Low et a l. 1989) and form vertical “ch imneys” in the inter - 
stellar medium (e.g. the Stockert ch imney. iMiiller et~aL 1 19871) . 
A general scenario was proposed by Norma n & Ikeuchi (1989) 
that chimneys play an important role in the interaction between 
the Galactic disk and halo by carrying up mass, energy, momen¬ 
tum, and magnetic flux. Therefore it is of interest to understand 
the properties of the superbubbles and their evolution processes. 
However, so far only a few superbubbles and chimneys have 
been resolved well and the magnetic field strengths measured 
(e.g. Vallee l 1993). 


The W4 superbubble in the Perseus arm is one of the few su¬ 
perbubbles that extend over several degrees and that have been 
also prev iously resolved we l l in ob servations. It was first iden¬ 
tified by iNormandeau et al.l (Il996l) from HI observations with 
the synthesis telescope of th e Dominion Radio Astrophysical 
Observatory (DRAO) ( Landecker et a lj 2000 ). A cone-shaped 
cavity was revealed to open upwards from its powering source, 
the young open star cluster OC1 352, which includes nine O- 
type stars. The lower part of the cavity is bounded by the well- 
known Galactic HII region W4. Evidence was found for an 
energetic outflow towards high positive latitudes, because HI 
streams were seen to point away from OC1 352. Because of the 
open conical shape viewed in the neutral gas, the entire struc¬ 
ture was first named “W4 chimney”. Lollow-up H a observations 
([Dennison et al.l!T997l) showed that the ionized gas is detected in 
the periphery of the HI cavity and revealed an enclosed super¬ 
bubble rather than a chimney, which is sealed at about b = 7°. As 
shown bv [West et al.l ( 2007 ). the appearance of this huge struc¬ 
ture seen from the 1.4 GHz radio continuum emission also re¬ 
sembles a closed bubble. We therefore refer to it as “W4 super¬ 
bubble” in this work. iBasu et af. ( 199 9) explained the morpho¬ 
logical difference between a chimney and a superbubble by the 
penetration of UV photons to the higher latitude regions. Based 
on the study of infrared and radio continuum data. lTerebev et al.l 
( 2003 ) find a large leakage of UV photons from OC1 352. This 
leakage was indirectly confirmed by the discovery of an even 
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Table 1 . Observational parameters 


Data 

Urumqi 46 cm 

Effelsberg 411 cm 

EMLS 421 cm 

Frequency [MHz] 

4800/4963 

2639 

1400 

Bandwidth [MHz] 

600/295 

80 

20 

HPBW['] 

9.5 

4.3 

9.35 

Number of coverages 

6(B)+4(L) 

6(B)+6(L) 

1(B)+1(L) 

Integration time [s] 

7.5 

6 

2 

rms (I/U, 0[mK T b \ 

0.7/0.4 

2.0/1.5 

15/8 

Calibrator 

3C286 

3C286 

3C286 

Flux density [Jy] 

7.5 

11.5 

14.4 

Polarization percentage 

11.3% 

9.9% 

9.3% 

Polarization angle [°] 

33 

33 

32 


larger high-latitude H a loop on top of the W4 superbubble ex¬ 
tending up to b = 30° (iRevnolds et al.ll200Th . 


The dist ance to the W4 superbubbl e was taken as 2 2 kpc i n 
the work of [N ormand eau et all (Il996l). IDennison et al.l (Il997l) . 
IRevnolds et ah ( 200ll). andlTerebev et all (20031). whil e 2.35 kpc 
was used bv lBasu et al.l (1 19991) and I West et aLN2007l) . Through 
measuring the trigonometric parallax of the methanol maser, 
IXu et al. ( 200 6) got 1.95 kpc for the distance of the W3 OH 
region. Considering a probable association between the HII re¬ 
gions W3 and W4 that gives another constraint on the distance 
of the W4 superbubble, we adopt the median value of 2.2 kpc in 
the following analysis. 


The age of the W4 jmp erbubble has not been determined 
with high precision. IDennison et a l. C1997) estimated an age be¬ 
tween 6.4 and 9.6 Myr, while iBasu et al.l (1 19991) derived a lower 
value of about 2.5 Myr. According to the formation of the large 
high-latitude H a loop found on top of the W4 superbubble, 
Reynolds et al.] (120011) argue that the formation of the W4 su¬ 
perbubble is not a single event, but rather caused by sequential 
star formation. Age estimates, however, are inevitably affected 
by the presence and strength of the magnetic fields. Numerical 
simulations have indicated that the magnetic field has a signif¬ 
icant impact on the evolu t ion of superbub bles (e.g. iTomisakal 
Il990t iFerriere et al.l Il99ll: ITomisakal 1 19981) . The expansion of 
superbubbles differs along and perpendicular to the _field lines 
and results in an elongated shape (Ide A villez & Brei tschwerdtl 
l2005l:lstn et al.ll2009h . lKomfienovic et al.l (Il999l) have suggested 
that magnetic fields of a few /iG must exist in the shells of the 
W4 superbubble to maintain its highly collimated shape and 
prevent its fragmentation from Rayleigh-Taylor instability. The 
strength of the line-of-sight component of the magnetic fields 
within the shells of the W4 superbubble was first estimated by 
[West et all ( 2007 ) based on the polarization data collected by 
the DRAO synthesis telescope. They determined a field strength 
between 3 /i G and 5 juG by analysing lines through the W4 sup- 
perbubble shell separated by 5'. No short-spacing information 
was included in their data, which is, however, required to cor¬ 
rectly interpret polarization structures resulting from Faraday 
rotation in the inter stellar medium (e.g. iReichl 12006 : ISun et al.l 
120071: lLandecker et"ai1l2010l) . We discuss the influence of polar¬ 
ized large-scale emission on this result in Sect. 4.2. 


This work presents a study of the radio spectrum and the 
magnetic fields within the shells of the W4 superbubble based on 
new sensitive multi-frequency radio continuum and polarization 
observations including zero spacings. We describe the data sets 
we used in Sect. 2 and present all observational results and their 
analysis in Sect. 3. In Sect. 4, we use a Faraday screen model to 
derive the magnetic field strengths in the W4 superbubble. We 
discuss and summarize the results in Sect. 5. 


2. Data 

New sensitive radio continuum and linear polarization observa¬ 
tions of the W4 superbubble have been made at 46 cm with 
the Urumqi 25-m radio telescope of Xinjiang Astronomical 
Observatories, Chinese Academy of Sciences, and at 411 
cm with the Effelsberg 100-m radio telescope of the Max- 
Planck-Institut fur Radioastronomie, Germany. Total intensity 
and polarization data at 421 cm were taken from an un¬ 
published sectioi^_o£_the Effel s berg Medium Latitude S urvey 
(EMLS) (lUvamker et al.l I1998L 119991: iReich et al.l l2004h. The 
Stockert421 cm total intensity northern sky survey (Reich 1982; 
IReich & Reichll 19861) and the DRA O 26-m single-dish polariza¬ 
tion survey (fWollebe n et al.l 120061) provided the missing large- 
scale intensities and the absolute zero levels. We also added 
high-resolution 421 cm dat a observed w ith the DRAO s ynthe- 
sis telescope by IWest et all (120071) and lLandecker et al.l (1201 Oh 
to the large-scale corrected EMLS maps. 

Technical details of each observation, such as the cen¬ 
tral frequency, angular resolution, and sen sitivity, are listed in 
Table d] The bandwidth depolarization (e.g. lCrawford & Tiffanvl 
120071 ) for the single-channel observations by the Urumqi and 
Effelsberg telescopes was found to be negligible (<1%) for all 
the three wavelengths, even for a rotation measure (RM) as large 
as 300 rad m -2 . 


2.1. Urumqi A6 cm data 

The 46 cm total intensity and linear polarization data were ob¬ 
tained between 2007 and 2009 during the conduction of the 
Sino-German 46 cm polarization s urvey of the Ga l actic plane 
( Sun et al] 120071: iGao et al.1 120101: ISun et al] 120111: IXiao et all 
l201ll) . using the same observational and reduction methods. The 
46 cm observations had two modes - the broad-band mode with 
a central frequency of 4 800 MHz and a bandwidth of 600 MHz 
- and a narrow-band mode with a central frequency shifted to 
4 963 MHz and a reduced bandwidth of 295 MHz. Raster scans 
with a speed of 47min were made along the Galactic longi¬ 
tude (L) and latitude (B) directions. A single observation mostly 
lasted for about two hours, so that only a portion of the large 
target field was covered. Elevations were always high enough 
that ground radiation fluctuations can be largely avoided. Ten 
full coverages were processed and then calibrated by 3C286 and 
3C295, the main polarized and un-polarized calibrators. Finally, 
individual maps were combined following the procedures de¬ 
scribed bv IGao et al . ( 12010 ). The final image of the W4 super¬ 
bubble (see Fig. [I]) is centred at l - 134?5, b = 5°6, covering an 
area of 5?3 x 5?7. 
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6cm TP (color) + 6cm TP contour 


6cm PI (color) + 6cm TP contour 


6cm PA (color) + 6cm PI contour 
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Fig. 1. Total intensity /, polarization intensity PI, and polarization angle PA images of the W4 superbubble at 46 cm ( upper panels 
a, b, c). All cm ( middle panels d, e, f), and A21 cm ( bottom panels g, h, i). The angular resolutions are 9'5 for A6 cm images in 
panels a and b, 4'3 for Al 1 cm images in panels d and e , and 9'35 for All cm images in panels g, h, and i. The PA images at A6 cm 
(panel c ) and Al 1 cm ( panel j) have angular resolutions of IT and 6', respectively. The total intensity contours for I and PI images 
run from 1.8 mK T t, in steps of 2.4 mK T t, for A6 cm, from 10.0 mK T^ in steps of 12.0 mK T& for Al 1 cm (after subtracting strong 
point-like sources), and from 2000 mK T^ in steps of 50 mK T^ for Al 1 cm. The PI contours on the PA images start from 5.0 mK 
Tfc in steps of 1.0 mK T& for A6 cm, and from 50.0 mK T& in steps of 8.0 mK T^ for 411 cm. PI contours are not overlaid on the 
421 cm PA image, because no clear correlations can be found. The rectangles in panel a are the regions for TT-plots study, while 
the rectangles in panels b, c, and e, /indicate the regions used for the Faraday screen model fitting. 


2.2. Effelsberg All cm data 

The 411 cm observations were done in the summer of 2008 with 
the Effelsberg 100-m telescope. The receiving system was de¬ 
scribed by lUvamker (2004 ), but was upgraded in 2005 with 
lower noise amplifiers and an eight-channel polarimeter. Here 
we made use of its broad-band channel with 80 MHz bandwidth. 
Data processing and calibration follow the same procedures as 
for the 46 cm observations. 


2.3. Combined A21 cm data 

The 421 cm data were extracted from the EMLS (Reich et al., in 
prep.) with a central frequency of 1 400 MHz and a bandwidth 
of 20 MHz. The missing large-sc ale structures were added by 
the S tockert 421 cm survey data dReichl[l982t iReich & Reichl 
1 19861) for total intensi ties, and the DRAO 26-m single-dish data 
( Wolleben et al.ll2006l) for polarization intensities. 

We also added the available 421 cm total intensity data from 
the Canadian Galactic plane survey dTavlor et al.1120031) . the po- 
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larization data from lLandecker et al.l (l201Qh . and the synthesis 
data from lWest et al.l (120071) . convolved to a 1'5 circular beam to 
the zero-level restored EMLS maps. A detailed description of the 
merging process of polarization data from the DRAO synthesis 
telescope, the Effelsberg 100-m and t he DRAO 26-m telescopes 
was given bv lLandecker et al' ( 2010 ). whose procedure we fol¬ 
lowed. The resulting T21 cm map does not fully cover the area 
of the A6 cm, All cm, and EMLS A21 cm maps. 


2.4. Absolute zero-level restoration for the A 6 cm and A 11 
cm polarization data 

Interferometric data miss short-spacings and single-dish maps 
are generally set on arbitrary zero levels. As emphasized by 
iReichl (120061) . missing large-scale structure means that diffuse 
polarized emission originating in the magnetized interstellar 
medium cannot be properly interpreted. 

The observed polarization U and Q maps at A6 cm and All 
cm were set to zero at their boundaries and thus miss polarized 
emission from larger-scale components. The zero-level problem 
was sol ved for the A6 cm polarization d ata of the Urumqi sur¬ 
vey (e.g. lSun et al.ll2007l:lGao et al.ll2010h by adding the missing 
large-scale component by extra polation of the WMAP K-band 
(22.8 GHz) polarization data dPage et al.l 120071: iHinshaw et al.l 
2009), which are on absolute levels. This is based on the as¬ 
sumption that Faraday rotation of the Galactic diffuse interstel¬ 
lar medium is negligible for frequencies as high as 4.8 GHz, 
which is valid except for the inner Galaxy dSun et al.ll201ll) . For 
the outer region within a few degrees away from the Galactic 
plane, such as the W4 superbubble, Faraday rotation effects can 
certainly be neglected between A6 cm and the WMAP K-band. 

For the W4 superbubble, we followed the same polariza¬ 
tion zero-level restoration scheme as described by iGao et al] 
( 20101). but using th e most recent WMAP K-band 9-yr data 
( Bennett et al.l 120131) . First, the Urumqi A6 cm polarization U 
and Q maps and the corresponding WMAP K-band polarization 
maps were convolved to a common angular resolution of 2°. The 
convolved WMAP U and Q maps needed to be scaled to A6 cm 
by a factor of where p is the brightness temperature 

spectral index (T& ~ V s ) for diffuse polarized emission. Spectral 
index maps for polarization intensities at 75' resolution were de¬ 
rived from the combined A21 cm (Effelsberg 100-m and DRAO 
26-m) and the WMAP K-band data, and also for control for the 
WMAP K- and Ka-band (33 GHz) data for the entire W4 super¬ 
bubble area. In both cases, the mean of was close to -3, i.e. 
for A21 cm versus K-band, p = -2.97+0.04, and for K-band ver¬ 
sus Ka-band, p = -2.98+0.31, indicating that depolarization at 
A21 cm is not important for large-scale emission. Finally, the fi- 
scaled WMAP U and Q data at 2° angular resolution were com¬ 
pared with the convolved A6 cm U and Q data, respectively. The 
differences were taken as the missing large-scale components 
and added to the observed A6 cm data at their original angular 
resolution. The A6 cm polarization intensity (PI) was then cal¬ 
culated from the restored U and Q as PI = ^/U 1 2 + Q 2 , and the 
polarization angle (PA) was obtained by PA = ^atan {U IQ). 

The Faraday rotation effect at Til cm is about 3.4 times 
that of A6 cm (A 2 dependence), so a successful zero-level re¬ 
covery using high-frequency WMAP K-band data requires a 
careful check on the rotation measure (RM) in the area of in¬ 
terest. For the W4 superbubble area, we calculated a RM map 
between T21 cm (Effelsberg + DRAO) and T1.3 cm (WMAP K- 
band) from the polarization angle data: RM = (PA i.3 Cm -PA2i C m) 
/ (0.013 2 - 0.214 2 ) at an angular resolution of 75'. The RMs 



137 ° 136 ° 135 ° 134 ° 133 ° 132 ° 


Galactic Longitude 

Fig. 2. Combined T21 cm polarized emission of the W4 super¬ 
bubble with data from the DRAO interferometer, the Effelsberg 
100-m, and the DRAO 26-m single-dish telescopes at an angular 
resolution of 3'. Overlaid T21 cm contours are from the source 
subtracted and spatially filtered total intensity map (see text). 
Contour lines run from 2000 mK T& in steps of 50 mK T&. 


are found to be small in general and the mean is 0.74 +1.01 
rad m -2 , which causes only small deviations in PA of -1° to 3°3 
between the Effelsberg Til cm and the WMAP T1.3 cm data. 
Therefore, we restored the missing large-scale polarized emis¬ 
sion for the Effelsberg Til cm data by the same procedures as 
for the Urumqi T6 cm data. 


3. Results 

We present the total intensity I, the polarization intensity PI of 
the W4 superbubble at T6 cm, Til cm, and T21 cm at angular 
resolutions of 9'5, 4'3, and 9'35 in Fig. [TJ where the T21 cm 
images are from the combined EMLS and DRAO 26-m single 
dish map. The angular resolutions of the T6 cm and Til cm PA 
images are convolved to 12' and 6', respectively, to increase the 
signal-to-noise ratio. The T21 cm PA image is shown at its orig¬ 
inal resolution of 9'35. A higher angular resolution T21 cm PI 
map combined from the DRAO synthesis telescope (IWest et ali 
120071) . the EMLS, and the DRAO 26-m telescope convolved to a 
3' circular beam is displayed in Fig. [2 This map is overlaid with 
contours of a source-subtracted and unsharp-masked T21 cm to¬ 
tal intensity map, which was combined from DRAO synthesis 
telescope, EMLS, and Stockert survey data. The singled-out PI 
images for the westerrQ shell of the W4 superbubble at T6 cm, 
Til cm, and T21 cm are shown in Fig. [3j For all the images, 
the T6 cm and Til cm total intensity images are the observed 
ones on relative zero levels, while all the others are restored to 
absolute zero levels. 


1 The rotation angle between the Galactic plane and the Celestial 

equator is small (~ 25°) in the W4 superbubble area, we thus refer to 
the larger longitudes as “east”, and the smaller longitudes as “west” in 
this work. 
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Fig. 3. Depolarization along the western shell of the W4 superbubble at 46 cm (panel a ), where total intensity contours start at 
4.9 mK Tfc and run in steps of 3.0 mK T&, and at 411 cm (panel b) where total intensity contours start from 9.0 mK T& in steps of 
13.5 mK T&, and at 421 cm (panel c) where the total intensity (DRAO interferometer + Effelsberg + DRAO 26-m) contours run 
from 2000 + (n-1) x 50 mK T& (n = 1, 2, 3, 4), and 2500 + (n-5) x 800 mK (n = 5, 6, 7, 8). The angular resolutions for the 46 
cm, 411 cm, and 421 cm images are 9'5, 4'3, and 3', respectively. 








Fig. 4. TT-plot results for the western shell (upper panels) and the eastern shell (lower panels ) of the W4 superbubble between 46 
cm and 411 cm (left panels) and for the four regions (A, B, C, and D) in the shells (see Fig. |T]) between 46 cm and 421 cm. 


3.1. Total intensity emission and spectral information 

In all total intensity images, the W4 superbubble shows up by 
two bright ridges emerging from the HII region complex W4 
at about i ~ 133° and i ~ 136°, extending towards high lati¬ 
tudes. These two ridges form a loop structure, since they close at 
about b - 6°. From the 46 cm and 411 cm images, a faint high- 


latitude extension of the western shell can be traced to more than 
b = 7°. Another partial shell-like structure can be distinguished 
above the major eastern shell of the W4 superbubble and merges 
with the western extension at about b - 8°. According to the ap¬ 
parent size determined by the maximum lengths on the Galactic 
longitude and latitude dimensions, AfxAb = 3° 8 x T.l (A b is 
calculated with respect to the position of OC1 352 at b - 029), 
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the physical size of the W4 superbubble is about 150 pc in width 
and 270 pc in height for a distance of 2.2 kpc. These results are 
consistent with previous estimates (e.g. [Dennison et all 19971) . 

Based on the tight correlation between radio continuum 
and H a emission and the flat radio spectrum derived between 
408 MHz and 1420 MHz, [West et al.l (120071) conclude that the 
radio continuum emission from the shells of the W4 super¬ 
bubble is optically thin thermal emission. We confirmed this 
by deriving the spectral index from the present data using TT- 
plots dTurtle et al.Hl 962l) betwe en 1.4 GHz and 4. 8 GHz. After 
subtracting the NVSS sources (ICondon et al.lll998l) . the bright¬ 
ness temperature spectral index between A6 cm and All cm 
was found to be p = -1.98 ± 0.16 for the western shell and 
P = -2.04 ± 0.26 for the eastern shell (Fig. ffl left panels ). 
Between A6 cm and A21 cm, we obtained ft = -2.12 + 0.39 for 
region ‘A’, p- -2.19 + 0.38 for region ‘B’, and/3 = -2.16 + 0.33 
and p = -2.16+0.10 for the regions ‘C’ and ‘D’, respectively for 
the marked regions in Fig. |T] These results agree with optical- 
thin thermal emission, which we use for the passive Faraday 
screen model to derive the magnetic fields within the shells of 
the W4 superbubble in Sect. 4.2. 

3.2. Polarization 

Unlike the resemblance in the total intensity images shown for 
the three different wavelengths, polarization behaves much dif¬ 
ferently. At A6 cm, the most pronounced feature related to the 
W4 superbubble is the strong depolarization that nicely follows 
the western shell and the PA deviation across the shell. No other 
prominent and closely related PI and PA features are seen along 
any other shell structures. 

The Al 1 cm polarization image discloses more structural de¬ 
tails owing to its higher angular resolution. The polarization in¬ 
tensities and angles vary strongly in the western shell area of 
the W4 superbubble (see Fig. [3]). Other than at A6 cm, depolar¬ 
ization is also seen for some sections of the eastern shell and the 
high-latitude filament centred approximately at l - 136 °,b = 7°. 
This, in conjunction with the A6 cm PI image, implies that RM 
varies in values and regularity over the W4 superbubble and is 
smaller in the eastern shell than in the western shell. Polarization 
structures are also visible in the bubble’s interior. A notable fea¬ 
ture is the arc-like structure running from longitude 133° to 134° 
and latitude 6?0 to 4?5, which is also visible but fainter at A6 cm. 

At A21 cm, no strong structural coincidence can be iden¬ 
tified between total intensity and polarization at first glance. 
[West et afl ( 2007 ) saw the western and eastern shells of the W4 
superbubble in depolarization in their A21 cm map using DRAO 
synthesis data alone. When averaging their data for a 10' wide 
latitude stripe centred at b = 4°. 5, the polarization intensity for 
the eastern shell drops from a surrounding level of about 65 mK 
T \y to 40 mK Tb, and for the western shell it decreases from 
50 mK Tb to about 40 mK T^. When combining the DRAO 
synthesis telescope data with the Effelsberg and DRAO single¬ 
dish data (Fig.0, we found that this depolarization is much less 
pronounced than the pure DRAO synthesis telescope data. The 
reason is the much higher absolute polarization level of about 
500 mK T b for the eastern and 350 mK T^ for the western shell. 
Fluctuations of the Galactic polarized emission seem to mask the 
depolarization along the W4 superbubble shells. For a zoomed 
view of the western shell at 3' resolution (Fig. [3]), there is, how¬ 
ever, a morphological coincidence between total and polariza¬ 
tion intensity, indicating that related depolarization exists. 

On larger scales, an inclined broad stripe of strong polar¬ 
ized emission runs across the image from the lower left (south¬ 


east) to the upper right (north-west) corner, while depolarization 
patches are seen in the orthogonal areas of the image. Both the 
Effelsberg and the DRAO single-dish data independently show 
this feature at A21 cm. This general distribution is also seen 
in the A6 cm and Til cm polarization data that are corrected 
by large-scale emission based on the WMAP K-band data (see 

Fig.ffl). 

4. The magnetic field strength in the W4 
superbubble shell 

Regular changes in sign or values of RMs of extragalactic 
sources across a finite area often indicate the existence of 
a Galactic Faraday screen. Sufficient RM sources in a sky 
area can be used to determine the line-of-sight component 
of the magnetic fields w ithin such foreground screens (e.g. 
lHarvev-Smith et all 120111) . There are 40 sources with deter¬ 
mined RMs within the W4 superbubble area in the catalogue 
of IXu & Hanl (f201 4 ]fl s upplemented by the NVSS RM cata¬ 
logue ( Tavlo r et al] 120091) . Most of these RMs show negative 
signs, resulting from the orientation of the Galactic magnetic 
field running clockwise (viewed from the north Galactic pole) 
in this part of the Perseus arm (lHan et alll2006l) . Five RMs are 
positive, which might indicate small-scale variations or source- 
intrinsic properties. However, the present RM number density is 
too sparse to estimate the magnetic field for shell regions of the 
W4 superbubble. 

4 . 1. The Faraday screen model 

A passive Faraday screen is a magneto-ionized bubble, fila¬ 
ment, or sheet that does not emit polarized emission, but mod¬ 
ulates background polarized emission passing through it via the 
Faraday effect. A successful model for deriving the physical pa¬ 
rameters of these Faraday screens was developed by [Sun et ali 
( 2007 ) for screens detected in the Urumqi A6 cm survey. It was 
also applied in several follow-up studies (e.g. iGao et al.ll2010l: 
IXiao et ap201lh . In its simplest form, the model takes the obser¬ 
vational fact that the foreground and background Galactic mag¬ 
netic field of the screens have the same orientation and PA ~ 0°, 
i.e. parallel to the Galactic plane at high enough frequencies. We 
have verified in Sect. 3.1 that the W4 superbubble exclusively 
emits optically thin thermal emission. We now check whether 
foreground and background PA s around the W4 superbubble 
area meet the condition of PA ~ 0°. Starlight polarization ob¬ 
servations of the W4 superbubble area (tHeilesll2QQQh support the 
idea that the Galactic magnetic fields are aligned to the Galactic 
plane (PA = -1° + 7°) for a distance range from 30 pc to 
2.85 kpc. This is further supported by the low-frequency po¬ 
larization measurements for the large excessive “Fan” region 
dSpoelstrall 19841) . covering the W4 superbubble area, and also by 
the WMAP K- and Ka-band polarization data, where Galactic 
RMs in this part of the Gala xy are small, and the large-scale 
Galactic magnetic fields (IHan et alJ120061) will only marginally 
change PA s from their intrinsic values. In addition, we also 
carried out A6 cm s imulations in the dir ection of the W4 su- 
perbubble using the ISun et al.l (120081) and ISun & Reichl (1201 Ol) 
Galactic 3D-emission model, which is based on all-sky surveys, 
RMs, and thermal emission data and also includes realistic mag¬ 
netic field turbulence (ISun & Reichl2009i) . The simulation result 
is shown in Fig. 0 Foreground PA of the W4 superbubble at 
2.2 kpc is close to zero and differs by just about 1° from the 

2 http://zmtt.bao.ac.cn/RM/index.html 
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Fig. 5. Simulations in the direction l - 135 °,b = 6° showing 
accumulated and spatial averaged RM (including one cr errors), 
PA, PI, and I values at A6 cm as a function of distance from the 
Sun. 


PA of the background emission at distances of 8 kpc or more. 
We repeated the simulation for various directions within the W4 
superbubble area at A6 cm, and obtained very similar results. 
Simulations at All cm and 421 cm show very small PA dif¬ 
ferences as well. We concluded that the Galactic 3D-emission 
simulations also support that the assumptions required for the 
passive Faraday screen model are satisfied. 

The passive Faraday screen model uses three parameters: 
a depolarization factor /, which describes the reduction of the 
background polarization by the Faraday screen; the Faraday 
screen foreground polarized emission fraction c with respect 
to the observed polarized emission at an off-position away 
from the screen; and the polarization rotation angle within the 
screen if/ s , which is the parameter for calculating RM (RM = 
if/si A 2 ), and the line-of-sight component of the magnetic field 
(see Eq. 6). These three parameters fit two observables: (1) the 
ratio PI on IPI o& and (2) the angle difference PA on - PA 0 s, where 
“on” and “off” denote that the line of sight passes through (on) 
or not (off) the screen. The Faraday screen model equation is 
given in Eq.fllfsee lSun et al1l2007[ for a detailed derivation): 


f PI ™ 


Ploff 

PAon 



(1 - c ) 2 + c 2 + 2fc(l - c) cos 2 if/ s , 


- PA 0 ff = - arctan 


/(I ~ c) sin 2^ \ 

c +/(1 - c) cos 2 ifr s ) 


( 1 ) 


In the following we allow / to be in the range [0.00, 1.00] 
and c between 0.6 to 1.0, when taking the simulation result 
(Fig. 0 into account. These parameters vary in steps of 0.01. 
if/ s is searched for in the range [-90°, 90°] in steps of 1°. 


4.2. Analysis of the western shell 

The clearest Faraday effect for the W4 superbubble is the depo¬ 
larized western shell seen at A6 cm and All cm. With the mea¬ 
surements at these two wavelengths, we estimated the line-of- 
sight component of the magnetic field Bp based on the Faraday 
screen model described above for a section as indicated by the 
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Fig. 6. Upper panel : profiles of the A6 cm I (relative zero level), 
PI (absolute zero level), and PA (absolute zero level) averaged 
perpendicular to the western shell of the W4 superbubble within 
the area marked by the large white rectangular region in Fig. Q] 
The pixel size is 3'. The vertical dashed lines indicate the “on - 
off” boundaries of the Faraday screen. Lower panel : the same as 
the upper panel , but for the Effelsberg All cm data. 


large white rectangle in Fig. Q] Averaged A6 cm and All cm I, 
PI, and PA profiles perpendicular to the western shell were cal¬ 
culated and shown in Fig. [6] The shell direction is inclined, and 
the averaged stripes or columns run parallel to it. We labelled the 
result by the column number from east to west across the west¬ 
ern shell (see Fig. [6]). The typical character of a passive Faraday 
screen is seen at A6 cm; i.e., an apparent PI depression is identi¬ 
fied with its minimum close to the peak position of the total in¬ 
tensity I. The very small PA difference between “on” and “off” 
positions implies a large foreground fraction c of the polarized 
emission. 

To calculate the two observables, PI on /PIoff and PA on ~ 
PA 0 ff, the “on” and “off” positions of the Faraday screen must 
be determined first. Combining the profiles of the A6 cm and All 
cm data, we set the boundary of “on” and “off” positions of the 
Faraday screen as marked in Fig. [6] The column size of 3' corre¬ 
sponds to 1.9 pc for a distance of 2.2 kpc of the W4 superbubble. 
The A6 cm PI Q ff and PA Q s were calculated as the average values 
of PI s and PAs on the “off” positions (Column 1-8 and Column 
17 - 25, see Fig.^upper panel). The central eight (9th to 16th, 
Fig.[6]and Table 2) columns represent the interior of the Faraday 
screen and were used for model fitting. 

We assumed lcr observational uncertainty for the A6 cm ob¬ 
servations and used Eq. 1 for the Faraday screen model fitting. 
A series of (/, c, i// s ) results were obtained, which are consis¬ 
tent with the A6 cm observables PI on /Ploff and PA on ~ PA 0 ff. 
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Table 2. Faraday screen parameters for the western shell of the W4 superbubble. Column numbers refer to Fig.|6] The best fit (/, c, i// s ) solution has the least deviations from the 
A6 cm and Al 1 cm observations. We tabulate n e (1 cr error = 5%) and the magnetic field component B// for filling factor f He = 0.4 and 1.0. We also list n e and B// in case the 12th 
column represents the longest path length (see text). 


Column No. 

9 

10 

11 

12 

13 

14 

15 

16 

0 l,b ) 

best fit (f.c.ffrj) 
f 

c 

(deg) 

RM (rad nr 2 ) 

132?86, 4?76 
(0.95, 0.85, -19°) 
0.77+0.15 
0.86+0.02 
-18+3 
-80+13 

132781, 4?78 
(0.78, 0.85, -25°) 
0.77+0.14 
0.87+0.02 
-20+3 
-89+13 

132?76, 4?79 
(0.95, 0.91, -59°) 
0.70+0.19 
0.90+0.02 
-64+4 
-286+18 

132771, 4781 
(0.68, 0.87, -64°) 
0.80+0.13 
0.88+0.02 
-64+3 
-286+13 

132?67, 4? 82 
(0.51,0.85, -69°) 
0.54+0.28 
0.87+0.02 
-66+6 
-295+27 

132?62, 4? 84 
(0.61, 0.88, -68°) 
0.49+0.28 
0.88+0.02 
-65+8 
-290+36 

132?57, 4?86 
(0.39, 0.86, -28°) 
0.26+0.19 
0.88+0.02 
-19+10 
-85+45 

132?52, 4? 87 
(0.59, 0.85, -13°) 
0.82+0.12 
0.86+0.01 
-16+2 
-71+9 

6cm I obs (mK T b ) 

6.2 

7.4 

8.7 

9.4 

8.6 

6.5 

3.9 

2.3 

EM (pc cm -6 ) 

47 

56 

66 

71 

65 

49 

30 

17 

path length (pc) 

25 

27 

31 

37 

61 

52 

40 

23 

© 

ll 

<4 

n e (cm 3 ) 

2.2 

2.3 

2.3 

2.2 

1.6 

1.5 

1.4 

1.4 

Bu (mG) 

-4.6+0.8 

-4.5+0.7 

-12.3+1.0 

-10.9+0.7 

-9.1+1.0 

-11.2+1.5 

-4.8+2.5 

-7.0+1.0 

fn e = 1.0 

n e (cm 3 ) 

1.4 

1.4 

1.5 

1.4 

1.0 

1.0 

0.9 

0.9 

Bu OuG) 

-2.9+0.5 

-2.8+0.4 

-7.8+0.6 

-6.9+0.5 

-5.8+0.6 

-7.1+1.0 

-3.0+1.6 

-4.4+0.6 

if the 12th column is taken as 
the longest path length 









o 

ll 

<4 

n e (cm 3 ) 

1.9 

1.9 

1.9 

1.6 

1.6 

1.5 

1.4 

1.4 

Bu OuG) 

-3.9+0.7 

-3.8+0.6 

-10.4+0.9 

-8.0+0.5 

-9.1+1.0 

-11.2+1.5 

-4.8+2.5 

-7.0+1.0 

fn e = 1.0 

n e (cm 3 ) 

1.2 

1.2 

1.2 

1.0 

1.0 

1.0 

0.9 

0.9 

Bu (JJ.G) 

-2.5+0.4 

-2.4+0.4 

-6.6+0.5 

-5.1+0.3 

-5.8+0.6 

-7.1 + 1.0 

-3.0+1.6 

-4.4+0.6 
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We then used these qualified (/, c, if/ s ) to predict All cm PI on 
(PI = ^jU 2 + Q 2 ) and PA on (PA = \ atari based on the fol¬ 
lowing equation (for details see lSun et al.ll2007l) : 

( U on =/^bgSin2^ , 

1 Son = PI fg +fPIbgCOs2ll/ s . 1 j 


Here, F/fgii cm — ^^offjlcm X C and FVbg.llcin — PI' off,llcm x 
(1 - c), where dll cm PI 0 s t n cm was calculated as PI 0 &, n C m = 

Woff,6cm ( 4800 MHz ^ = 65 - 8 mK T *’ with a s P ectral index p = -3 
as obtained in Sect. 2.4. This agrees with the measured values 
at the edges of the dl 1 cm PI profile in Fig. [6j The rotation an¬ 
gles at dll cm were calculated as if/ncm = 6cm ( qo^ ) 2 - We 
compared the predicted values with the dll cm observational 
results within lcr errors, which further filtered the (/, c, if/ s ) data 
to agree with both d6 cm and dll cm observations. From all the 
remaining results, we found the parameter c to be rather well 
constrained within 0.8 to 0.9 (see Table [2]) and largely indepen¬ 
dent of / variations. This is consistent with the dominant fore¬ 
ground contribution towards the W4 superbubble. The reduction 
of the polarized background emission is mathematically possible 
in various ways. Here, if/ s is the most interesting parameter and 
is also well constrained, which enables us to make a reasonable 
estimate of the magnetic field and assess its errors. In Table 2, 
we listed the average values of /, c, and if/ s for each column in 
the Faraday screen, including standard deviation. We searched 
for the best fit (/, c, i// s ) for each column by the least x 2 method. 
The difference between the average and the best fit value (/, 
c, if/ s ) is small, with the maximum of 1.9<x. In the following, 
we used the averaged i// s and its standard deviation to derive the 
magnetic field within the western shell of the W4 superbubble. 

The T21 cm observations cannot be used for the model fit¬ 
ting since depolarization occurs at T21 cm even by smaller RMs 
within the shell. However, it can be examined to verify the fit re¬ 
sult. Considering the small-scale RM structures as shown by the 
Al 1 cm data, total depolarization could be expected for the back¬ 
ground PI at A2 1 cm passing through the bubble’s shell. In fact, 
the interferometric T21 cm image bv [West et al.l ( 2007 ) shows 
the depolarization along the W4 superbubble shells, which is al¬ 
most masked in the single-dish polarization image as discussed 
in Sect. 3.2. In the [West et al . (200 7) interferometric map, the 
depolarization averaged for a 10' wide latitude stripe is about 
25 mK Tb along the eastern and about 10 mK T^ along the 
western shell. The single-dish map polarization level of about 
500 mK T^ in the eastern and about 350 mK T& in the western 
shell area imply an upper limit of about 0.97 for c in the west¬ 
ern shell and about 0.95 in the eastern shell. The large difference 
in the polarization emission levels in the interferometric and in 
the single-dish maps means that a quantitative analysis of the 
interferometric polarization data without large-scale correction 
cannot be made to estimate the field strength. We concluded that 
the A2 1 cm polarization data agree with the parameters of the 
Faraday screen model along the western shell based on the A6 
cm and Al 1 cm observations. 

For an estimate of Bn within the western shell, the thermal 
electron density n e , its filling factor f„ e , and the path length / 
within the bubble’s shell are needed. Emission measure (EM), 
which is the integral of f He and n e 2 over the path length dl can 
be calculated, because it is related to optically-thin thermal ra¬ 
dio continuum emission by the following equation (IWilson et al.l 
12013h : 


(|, = 8.235 x,0-(|r»», GHz 




EM 
pc cm" 


-)a(v,T), (3) 


where, 7 \ is the radio continuum brightness temperature mea¬ 
sured at frequency v, a(v,T) is a factor close to 1. We assumed T e 
= 8 000 K for the ionized Galactic interstellar medium in general. 
Based on the A6 cm observation, we deduced the EM s within the 
western shell of the W4 superbubble. However, from Fig. [6] we 
saw a large-scale emission component of about 1.9 mK T& at 
A6 cm besides thermal emission I from the western shell. After 
subtracting this large-scale component, EM s of the central eight 
columns within the screen area were calculated based on Eq. 0 
as listed in Table 2. 

To verify the EM results that we obtained from the radio con¬ 
tinuum data, H a measure ments were used to estimate the EMs 
indep endently following Haffn er et al.l (1 19981) and iFinkbeinerl 
(120031) by using 

EM = 2.75 T 9 ' 9 I Ha exp[2.44£(£ - V)]. (4) 

iDennison et al.l ([1997) measured H a emission for several po¬ 
sitions towards the W4 superbubble. At their position K (l ~ 
133?0, b ~ 52 4), they observed an intensity of 5.8 Rayleigh. 
[Reynolds et al. ( 200 1) measured 6.9+0.5 Rayleigh at the posi¬ 
tion £ = 13228, Z? - 5°. 1. Both directions are within the western 
shell region of the W4 superbubble. The extinction parameter 
E(B - V), which must be taken into account, is in general diffi¬ 
cult to obtain and often ignored for regions outside the Galactic 
plane. Two stars, Hilt 266 (l = 13222,Z? = 3°.l) and Hilt 338 
(l = 134?2, b — 3?0), were found near the W4 superbubble with 
distances of 2.5 kpc and 2.3 kpc, respectively. They have a mea¬ 
sured extinction value of E(B - V) = 0.66 (tFfiltneif19561) . Based 
on Eq.[4]and assuming T = 8 000 K, we derived the EMs for the 
two positions as 65 pc cm -6 and 77.7+5.6 pc cm -6 , which are 
consistent with those based on radio continuum data in Table 2. 

The path lengths for the central eight columns (9th -16th) 
across the western shell of the W4 superbubble are related to the 
geometry, where we assume a uniform shell. The largest RM cor¬ 
responds to the longest path length dValleel 19821) . which is corre¬ 
sponding to the 13th column in our case (see Fig.[6]and Table 2), 
and defines the inner radius of the shell. We also considered the 
case of a thermal shell, where the strongest total intensity emis¬ 
sion is at the 12th column, and included the results in Table 2. 
The angle between the central axis of the W4 superbubble and 
the 13th column is about 1°7, and the angle between the 13th 
column and the outer edge of the western shell is 1(T5. Then for 
a distance of 2.2 kpc, we calculated the inner and outer radii of 
the shell as 65 pc and 72 pc, respectively. With a column size of 
3', the path lengths of the central eight columns were deduced. 
The electron densities and the strength of the line-of-sight com¬ 
ponent of the mag netic field Bn can be obtai ned following the 
equations given bv lHarvev-Smith et al.l (2 01 ll) . where the filling 
factor f ne is the only unknown parameter: 



(5) 


B/i = 


RM 

0.81 y/EMy/fiJ 


yUG. 


( 6 ) 


Th e fillin g factor is often assumed to be 1.0, as in I West etakl 
(120071) . while lHarvev-Smith et al.l (1201 ll) used f He =0.1 for the 
sample of HII regions they studied. Since no well-determined 
value has been reported for superbubbles so far, we calculated 
n e and B // based on Eqs. [5] and 0 for f„ e = 0.1, 0.4, 0.7, and 
1.0. We used the standard error propagation to estimate the er¬ 
rors in n e and B//. The uncertainty of n e results from the error of 


9 























































Gao et al.: W4 superbubble 


00 



a 




Fig. 7. Left panel : electron density n e across the western shell of the W4 superbubble for filling factors f„ e =0.1 (green triangle), 0.4 
(blue square), 0.7 (cyan diamond), and 1.0 (magenta circle). Right panel : the strength of the line-of-sight component of the B-field 
across the western shell for different filling factors f Ue . 


EM , which stems from the A6 cm total intensity measurements. 
ISun et al . ( 2007 ) quote the maximum uncertainty of 10% for the 
A6 cm survey data. Besides the uncertainty of EM, the errors in 
Bn also depend on the uncertainties in RM. The standard devia¬ 
tions used for this purpose were listed in Table 2. 

We showed the n e and Bp profiles across the western shell of 
the W4 superbubble in Fig. [71 Owing to the 1/ ^Jfn e dependence, 
fn e = 0.1 results in high electron density and a large B field, while 
the values for f He = 0.4, 0.7, and 1.0 are much smaller. To dis¬ 
tinguish which case is more realistic, we made a simple esti¬ 
mate according to the magnetic flux conservation. For the W4 
superbubble, the swept-up magnetic field should be a factor of 
about 5 higher than the undisturbed Galactic B field when we 
consider a ring with an inner radius of 65 pc and jm outer ra¬ 
dius of 72 pc. Using RM data from pulsars, lHan et ali ( 2006 ) 
find that the strength of the large-scale regular radial B field of 
the Galaxy decreases from the inner Galaxy (4 pG) to the solar 
neighbourhood (2 pG). If a turbulent field of the same order is 
considered, the field strength in the W4 s uperb ubbl e area is about 
3 pG. The 3-D simulation of ISun et ahl (120081) and lSun & Reichl 
( 2009) where both regular and turbulent B fields were included 
gave a similar value of B tot of around 2.4 pG. Thus the enhanced 
line-of-sight B-field strength in the western shell should not ex¬ 
ceed 15 pG by much, otherwise an amplification mechanism is 
needed to explain a larger B field for f He = 0.1. If a B-field ampli¬ 
fication can be excluded, the filling factor f Ue in the western shell 
of the W4 superbubble is likely within the range of 0.4 to 1.0. We 
limit the following discussion to the case that the filling factor f He 
is 0.4, 0.7, and 1.0, where f He = 1.0 gives the lower limit of the 
electron density n e and the strength of the line-of-sight compo¬ 
nent of the magnetic field B / /. When considering that the line- 
of-sight passes solely through the western shell (13th to 16th 
column), we obtained the average B// to be -5.0(±10%) pG (J He 
= 1.0) to -8.0(±10%) pG (f ne = 0.4), and the electron density n e 
is between 1.0(±5%) cm -3 (f„ e = 1.0) and 1.5(±5%) cm -3 (f„ e = 
0.4). 

[West et all ( 2007 ) used a different analysis method to de¬ 
rive the line-of-sight component of the magnetic field strength 
in the shell of the W4 supperbubble. Based on the observed de¬ 
polarization amount, they estimated the PA difference from ad¬ 
jacent lines passing through the shell. From the Faraday screen 
model, we found the RMs within the shell for the 13th to 16th 


column between -70 rad m -2 and -300 rad m ~ 2 , w hich causes 
one to four full rotations of PAs at A21 cm. IWest et all (120071) 
quote the A21 cm differential angle of 60° and calculated a line- 
of-sight B-field strength between 3.4 pG and 9.1 pG from the 
different path lengths through the shell. The polarization in¬ 
tensity, polarization angles, and their differences are, however, 
influenced by the large-scale polarized emission, which is in¬ 
completely observed by the DR AO interferometer. I Sun et ali 
( 2007 ) quote the observed Stokes parameter U , Q from a pas¬ 
sive Faraday screen as U = PIf g sin2if/f g +fPIb g sm(2iffb g + 2i// s ), 
Q = PIf g cos2ij/f g + fPIb g cos(2if/b g + 2if/ s ), where fg and bg 
denote foreground and background components of PI and PA. 
Here, / is the depolarization factor and i// s is the Faraday ro¬ 
tation angle within the screen. For the W4 superbubble, the 
foreground and background PA angles are close to zero (see 
Sect. 4.1.), then U and Q were calculated as U = fPIb g sin2i// s 
and Q = PIf g + fPIb g cos2if/ s (Eq. 2). Therefore, the observed 
depolarization amount in PI and the PA differences for neigh¬ 
bouring lines through the shell depend not only on the differ¬ 
ence in xf/ s , which is needed for a B-field calculation, but also 
on PIf g and Ph g - Without taking this into account, the B-field 
result depends in some way on the large-scale polarization frac¬ 
tion included in the data. Only if the Faraday screen is very local 
and PIf g = 0 mK T^ does this method give the correct results. 
However, this is not the case for the W4 superbubble. 

4.3. The north-eastern extension 

For better understanding of the magnetic field configuration and 
the strength of the W4 superbubble, we searched for more parts 
of shells and arcs, where Faraday effects could be studied in 
terms of the described Faraday screen model. We note that it 
is difficult to find unconfused areas with good quality “on” and 
“off” data. However, a section of the north-eastern extension lo¬ 
cated above the major eastern shell is noted in the Til cm po¬ 
larization image (see Fig. [8]), which we found to be suitable for 
modelling. The A6 cm data show no evidence of depolarization 
in this region, indicating that the RM is smaller than in the west¬ 
ern part of the W4 superbubble. We again applied the Faraday 
screen model and obtained / ~ 0.87, c ~ 0.86 and a positive 
RM of about 55 rad m -2 in this area. The fitted parameters pre¬ 
dict PI on /PI Q tf = 0.97 and PA on - PA 0 s = 1?5 for the A6 cm, 
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Fig. 8. /Ill cm PI (left panel) and PA (right panel) of the north-eastern shell extension overlaid by total intensity contours. The 
contour levels and the angular resolution of the images are the same as in Fig. [U 


and PI on I PI off = 0-89 and PA on - PA 0 ff = -3°9 for the A2\ cm 
observations, which generally agrees with the non-detection of 
noticeable Faraday rotated structures in both maps. Following 
the procedures introduced above, we obtained a lower electron 
density of n e = 0.5 (f„ e = 1.0) - 0.9 cm -3 (f„ e = 0.4) with 1 cr error 
of 15%, and a weaker line-of-sight component of the magnetic 
field B// = 3.1 ( f ne = 1.0) - 4.9 pG (f He = 0.4) with ler error of 
30% within the upper north-eastern shell extension. 

5. Discussion and summary 

We have shown that strong regular magnetic fields exist in 
the W4 superbubble’s shells. A study of the process of cre¬ 
ating this magnetic field configuration is beyond the scope of 
this paper, while numerical sim ulations (e.g. lTomisakal [l998: 
Ide Avillez & Breitschwerdtll2005l) provide hints of an evolution¬ 
ary scenario. To explain the observations, it seems most likely 
that the originally plane-parallel magnetic fields were lifted up 
and compressed during the expansion of the W4 superbubble. 

The calculation of the total magnetic field strength in the 
W4 superbubble shell depends on its orientation relative to the 
line of sight. From the disco very of atomic hydrogen asso¬ 
ciated to the W4 superbubble, iNormandeau et al.l ( 1996 ) con¬ 
cluded that the W4 sup erbubble is inclined towards the observer. 
lLagrois & Joncasl ( 2009 ) studied this inclination in detail and 
found that the inclination angle is between 9° and 27°. In ad¬ 
dition, we measured an inclination angle of 18° with respect 
to the Galactic north pole for the western shell (see Fig. [TJ. 
If the magnetic field in the W4 superbubble is directed along 
its shell and if taking its geometry into account, the depro- 
jected regular B tot field in the western shell can be deduced as 

B tot = B/i ll + (— 1 0 ) 2 = (2.3 - 6.7) Bn, where 6 is the 

inclination angle between 9° and 27°. 

Thermal and magnetic pressure, besides gravity, ram pres¬ 
sure, and turbulence are key input parameters for superbubble 
evolution simulations. Thermal pressure is defined as P t her - 


2nokT e , where no = n e for total ionization, k is the Boltzmann 
constant, and T e is 8 000 K as assumed above. Magnetic pres¬ 
sure is calculated as P mag = B^ ot /%n. For each column within 
the western shell of the W4 superbubble shell (13th, 14th, and 
16th columns, the 15th column is excluded because of the poor 
constraints on the B-field), we calculated the thermal and mag¬ 
netic pressure based on the results shown in Fig. [7] and made a 
comparison in Table 3. For the inclination angle of 6 = 27°, the 
magnetic pressure P mag is generally several times of the thermal 
pressure Pther, and P mag becomes overwhelming when the incli¬ 
nation angle 6 approaches to 9°. 

For the north-eastern extension in the higher latitude region, 
P mag and P t her are comparable when 6 ~ 27°, regardless of f Ue 
= 0.4 or 1.0. P mag dominates P t h er if 0 approaches 9°, even for 
fn e = 1.0. A more precise inclination angle would help to better 
constrain the physical properties of the W4 superbubble. 

We found a positive RM on the eastern side and a nega¬ 
tive RM on the western side of the W4 superbubble. This is 
expected for the scenario described above, where the W4 su¬ 
perbubble expands and breaks out of the Galactic plane and lifts 
up the magnetic field, which runs clockwise_ (viewed from the 
north Galactic pole) in the Perseus arm (lHan et al1l2006 ). Then, 
for the eastern shell, the field lines will go up and for the western 
shell downwards. Because the W4 superbubble til ts towards us 
(INormandeau et al.lIl996h lLagrois & Joncasl 120091) . the line-of- 
sight component of the magnetic field points away in the west¬ 
ern shell of the W4 superbubble, resulting in a negative RM, and 
towards us in the eastern shell, where a positive RM is observed. 

Radio continuum and polarization observations at A6 cm, 
All cm, and A21 cm have been made to study the radio emission 
properties and to estimate the magnetic fields of the W4 super¬ 
bubble. With the flat radio continuum spectrum found between 
A6 cm, All cm, and A21 cm, we confirm the thermal origin of 
the radio continuum emission of the W4 superbubble. Polarized 
emission shows dramatically morphological differences from 
wavelength to wavelength. With the advantage of being less af¬ 
fected by Faraday rotation, the A6 cm and Til cm polarization 
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Table 3. Thermal pressure versus magnetic pressure in the western shell of the W4 superbubble for different filling factors f Ue and 
the inclination angles of 0 = 27°, 6 = 18°, and 6 = 9°. 


Column No. 

Pr/te/ 

(10- 12 dyn cm" 2 ) 

p mos , e = 27° 
(10- 12 dyn cm" 2 ) 

v mag , e = 18 ° 

(IQ" 12 dyn cm" 2 ) 

P mag. 9 = 9° 

(10 -12 dyn cm" 2 ) 

fn e = 0.4 

13th 

3.6+0.2 

17.5+3.7 

38.2+8.1 

150.0+31.9 

14th 

3.4+0.2 

26.3+7.1 

57.4+15.4 

225.5+60.5 

16th 

3.0+0.2 

10.3+2.9 

22.4+6.3 

88.1+24.6 

fn e = 0.7 

13th 

2.7+0.1 

10.0+2.1 

21.8+4.6 

85.7+18.2 

14th 

2.6+0.1 

15.0+4.0 

32.8+8.8 

128.9+34.6 

16th 

2.3+0.1 

5.9+1.6 

12.8+3.6 

50.3+14.1 

fn e = 1.0 

13th 

2.3+0.1 

7.0+1.5 

15.3+3.2 

60.0+12.7 

14th 

2.1+0.1 

10.5+2.8 

23.0+6.2 

90.2+24.2 

16th 

1.9+0.1 

4.1+1.1 

9.0+2.5 

35.2+9.8 


data were used for the estimates of the line-of-sight component 
of the magnetic field within the western shell of the W4 super¬ 
bubble by a passive Faraday screen model. Considering the ther¬ 
mal electron filling factor f He , the radio continuum observations 
and a simple geometric assumption result in an electron density 
of n e ~ 1.0 x 1/ ^]J^ e cm -3 , and the line-of-sight component of 

the magnetic field B// -5.0x1/ jjG (pointing away from 

us) for the western shell of the W4 superbubble, where the typi¬ 
cal error is about 5% for n e and 10% for B//. Based on a simple 
estimate, we found that L likely has a value great er than 0.1 
assumed for HII regions (lHarvev-Smith et all 1201 lh . Being re¬ 
lated to the inclination angle of the superbubble with respect to 
the plane of the sky, a total magnetic field is found above 12 /iG. 
This results in a magnetic pressure, which is one or two orders of 
magnitude higher than the thermal pressure in the western shell 
of the W4 superbubble. The All cm polarization data allow a 
model fit of the weaker RM in the high-latitude north-eastern 
shell, where we find a lower electron density of n e = 0.5 to 0.9 
cm -3 and B// of 3.1 - 4.9 ji G for filling factors f He of 1.0 or 0.4, 
respectively. This means that the magnetic and thermal pressure 
might be comparable in the upper parts of the superbubble. The 
RM sign reverses as expected for a scenario where the Galactic 
magnetic field is pushed out of the plane by the expanding W4 
superbubble, which is tilted towards us. These results are ex¬ 
pected to constrain magneto-hydrodynamical simulations of the 
W4 superbubble and superbubbles in general. 
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